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Abstract 

In this paper, we study inverse scattering of massless Dirac fields that propagate in the exterior 
region of a Reissner-Nordstrom black hole. Using a stationary approach we determine precisely the 
leading terms of the high-energy asymptotic expansion of the scattering matrix that, in turn, permit 
us to recover uniquely the mass of the black hole and its charge up to a sign. 

1 Introduction 

Black hole spacetimes are probably among the most fascinating objects whose existence is predicted 
by Einstein's general relativity theory. In the last twenty years, numerous mathematical studies have 
been achieved in order to better understand their properties. Propagation of fields (wave, Klein-Gordon, 
Dirac, Maxwell) in these peculiar geometries and direct scattering results have been extensively studied 
in [U [31 [71 [H Ell El EH1 H2] and used for instance to give rigorous interpretation of the Hawking effect 
[H [T8l [23] . Other original phenomena related to black hole spacetimes are for instance superradiance 
(the possibility of extracting energy from a rotating black hole) [3 [12l [13l [14] or causality violation (the 
existence of time machine inside rotating black hole that allows to travel through time) [6]. Despite the 
richness of these original phenomena, a striking feature of black holes spacetimes is their simplicity. Here 
we refer to the fact that, whatever might be the initial configuration of a collapsing stellar body, the 
resulting black hole spacetime can be eventually described by "at most" three parameters - its mass, 
its electric charge and its angular momentum - an important uniqueness result summarized by the well- 
known formula: Black holes have no hair. A natural question arises thus: assume that we are observers 
living in the exterior region of a black hole, at rest with respect to it and located far from it (such 
observers are called static at infinity and can be typically thought as a telescope on earth aimed at the 
black hole), can we in these conditions measure the defining parameters of the black hole by an inverse 
scattering experiment? 

The result contained in this paper is a first step in this direction. Here we focus our attention on 
Reissner-Nordstrom black holes, that is black holes that are spherically symmetric and electrically charged 
and thus described by only two parameters: their mass M and their charge Q. We consider massless 
Dirac fields which evolve in the exterior region of a Reissner-Nordstrom black hole and we use the direct 
scattering results, already obtained in [24j [22l [Jj , to define the wave operators and the corresponding- 
scattering operator S. The main result of this paper is then the following. Suppose that the scattering 
operator S is known to observers static at infinity (accessible in theory by physical experiments), then 
we show that such observers can recover uniquely the mass M of the black hole as well as its charge Q, 
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up to a sign. Note here that we cannot completely determine the charge Q of the black hole because 
of the simplicity of our model. The propagation of massless Dirac fields indeed is only influenced by 
the geometry {i.e. the metric) of the black hole which depends on Q 2 (cf. (|2.1[) ). The consideration of 
massive charged Dirac fields should solve this undetermination since an explicit term corresponding to 
the interaction between the charge of the black hole and that of the field appear then in the equation (see 
[7]). The choice of dealing with massless Dirac fields is nevertheless meaningful since it is the simplest 
model we could study (in particular no long-range terms appear in the equation) that already provides 
significant informations. Moreover, we expect to be able to extend this result to much more complicated 
models such as rotating black holes (or Kerr black holes) for which there already exists a direct scattering 
theory [19] . 

The strategy we adopt to prove our main result is based on a high-energy asymptotic expansion of the 
scattering operator S. Such a technique was introduced by Enss and Weder in [11] and used sucessfully to 
recover the potential of multidimensional Schrodinger operators (note that the case of multidimensional 
Dirac operators in flat spacetime was treated later by Jung in 21j). They showed that the first term of 
the high-energy asymptotics is exactly the Radon transform of the potential they are looking for. Since 
they work in dimension greater than two, this Radon transform can be inversed and the potential thus 
recovered. In our problem however, due to the spherical symmetry of the black hole, we are led to study 
a one dimensional Dirac equation and the Radon transform cannot be inversed in this case. Nevertheless 
the first term of the high-energy asymptotics leads to an integral that can fortunately be explicitely 
computed and gives exactly the radius of the event horizon, already a physically relevant information. 
In order to obtain the mass and the charge of the black hole, we must then calculate the second term 
of the asymptotics. We follow here the stationary technique introduced by one of us [25] which is close 
in spirit to the Isozaki-Kitada method used in long-range scattering theory [20]. The basic idea is to 
replace the wave operators (and thus the scattering operator) by explicit Fourier Integral Operators, 
called modifiers, from which we are able to compute the high-energy asymptotic expansion readily. The 
construction of these modifiers and the precise determination of their phase and amplitude is the crux 
of this paper. While this method was well-known for Schrodinger operators and applied successfully to 
various situations (see [1, 25, 26, 27 ), it has required some substantial modifications when applied to our 
model due to the specific properties of Dirac operators, essentially caused by the matrix- valued nature 
of the equation. A good starting point to deal with this difficulty was the paper of Gatel and Yafaev 
[15] where direct scattering theory of massive Dirac fields in flat spacetime was studied and modifiers (at 
fixed energy) were constructed. 

2 Reissner-Nordstrom black hole and Dirac equation 

In this section, we first briefly describe the geometry of Reissner-Nordstrom black holes. In particular 
we point out the main features and properties that are meaningful from scattering theory viewpoint. We 
then write down the partial differential equation that governs the evolution of massless Dirac fields in 
the exterior region of Reissner-Nordstrom black holes and recall some known results in direct scattering 
theory. 

2.1 Reissner-Nordstrom black holes 

In Schwarzschild coordinates a Reissner-Nordstrom black hole is described by a four dimensional smooth 
manifold 

M=R t xl+x Si, 



2 



equipped with the lorentzian metric 

g = F(r)dt 2 - (F(r)) _1 dr 2 -r 2 diu 2 , (2.1) 

where 

and gL> 2 = <i6> 2 + sin 2 ci^ 2 is the euclidean metric on the sphere S 2 . The quantities M > and Q e K 
appearing in (|2.2[) are interpreted as the mass and the electric charge of the black hole. 

The metric (|2.1j) has two types of singularities. Firstly, the point {r = 0} for which the function F is 
singular. This is a true singularity or curvature singularity 0- Secondly, the spheres whose radii are the 
roots of F (note that the coefficient of the metric g involving explodes in this case). The number of 
these roots depends on the respective values of the constants M and Q. In this paper we only consider 
the case M > \Q\ for which the function F has two zeros at the values r± = M± \J M 2 — Q 2 . The spheres 
{r = 7' + } and {r = r_} are called the exterior and interior event horizons of the Reissner-Nordstrom 
black hole. Both exterior and interior horizons are not true singularities in the sense given for {r = 0}, 
but coordinate singularities. It turns out that, using appropriate coordinate systems, these horizons can 
be understood as regular null hypersurfaces that can be crossed one way but would require speeds greater 
than that of light to be crossed the other way. Hence their name: event horizons (we refer to [30] for a 
general introduction to black hole spacetimes). 

Throughout this paper we shall nevertheless use the Schwarzschild coordinates as the coordinates of 
our analysis. The reasons are twofold. First in this coordinate system the coefficients of the metric (|2.ip 
don't depend on the variables t and u>, reflecting in fact the various symmetries of the black hole. We 
shall see in the next subsection 12.21 that the Dirac equation will take a very convenient form, namely 
an evolution equation (with respect to t) with time-independent Hamiltonian acting on the spacelike 
hypersurfacc £ = M+ x S 2 , a nice formalism when dealing with scattering theory. Second, it corresponds 
implicitely to the natural notion of observers static at infinity mentioned in the introduction. Such 
observers are located far away from the exterior horizon of the black hole and live on the integral curves 
of their velocity 4-vector 

U = — — — 

F(r) dt 

In particular, the time coordinate t is the proper time of such observers and thus corresponds to their 
true experience of time. We refer to [28 for a more complete discussion of this notion.. This choice of 
coordinates appears then quite natural with the idea of scattering experiments we have in mind. 

We point out now a remarkable property of the exterior horizon when described in Schwarzschild 
coordinates, property that has important consequences for scattering theory. From the point of view 
of observers static at infinity, the exterior horizon is perceived as an asymptotic region of spacetime. 
Precisely, it means that the exterior horizon is never reached in a finite time by incoming and outgoing 
null radial geodesies, i.e. the trajectories followed by light-rays aimed radially at the black hole or at 
infinity. These geodesies are given by the integral curves of F(r)~ x ^ ± Along the outgoing geodesies 
we can express the time t as a function of r by the formula 

t(r) = f F(r)- l dr, (2.3) 

where ro > r + is fixed. From (|2.2[) and (|2.3[) . we see immediately that t(r) — > — oo when r — > r+. An 
analogue formula holds for incoming geodesies. 



x It means that certain scalars obtained by contracting the Riemann tensor blow up when r — > 0. 
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In consequence we can restrict our attention to the exterior region {r > r + } of a Reissner- Nordstrom 
black hole and study the inverse problem for massless Dirac fields there. We just have to keep in mind 
that we won't need any boundary conditions on the exterior event horizon {r = r + } since this horizon is 
an asymptotic region spacetime. To make this point even clearer, let us introduce a new radial coordinate 
x, called the Regge- Wheeler coordinate, which has the property of straightening the null radial geodesies 
and will greatly simplify the later analysis. It is defined implicitely by the relation 



£-m (2-4) 



and explicitely by 



x = r H — log(r — r + ) H log(r — r_). (2.5) 

r+ — r_ r+ — r_ 

In the coordinate system (t,x, lj), the horizon {r = r+j is pushed away to {x = —00} and thanks to 
(|2.4[) . the metric takes the form 

g = F(r){dt 2 -dx 2 )-r 2 du?. (2.6) 

Observe that the incoming and outgoing null radial geodesies are now generated by the vector fields 
± -gz and take the simple form 

7 ± (t) = (t,xo±t,u Q ), <GM, (2.7) 

where (xo,u>q) elxS 2 are fixed. These are simply straight lines with velocity ±1 mimicking, at least 
in the t — x plane, the situation of a one-dimensional Minkowski spacetime. 

From now on we shall work on the manifold B = M. t x R x x S 2 equipped with the metric (12. 6|) . 

2.2 Dirac equation and direct scattering results 

Scattering theory for Dirac equations on the spacetime B has been the object of several papers [24l [22l [7] 
(chronological order) in the last fifteen years. We use the following convenient form for a massless Dirac 
equation obtained therein. Under Hamiltonian form the equation eventually reads 

id t 4> = H^j, (2.8) 

where -0 is a 2-components spinor belonging to the Hilbert space H = L 2 (M. x S 2 ; C 2 ) and the Hamiltonian 
H is given by 

H = T 1 D x + a(x)D S 2, (2.9) 
where a(x) — , D x — —id x , D32 denotes the Dirac operator on S 2 , i.e. 

D S , = - t T 2 (de + A - — n T 3 d v , (2.10) 
2 sin 

and r 1 , r 2 , r 3 appearing in (|2.9p and (|2.10[) are usual 2x2 Dirac matrices that satisfy the anticommutaton 
relations PP + T^T 1 = 2Id. 

We can simplify further the expression of the Hamiltonian by using the spherical symmetry of the 
equation. The operator D52 has compact resolvent and can be diagonalized into an infinite sum of 
matrix- valued multiplication operators. The eigenfunctions associated to D52 are a generalization of 
usual spherical harmonics called spin-weighted harmonics. We refer to I.M. Gel'Fand and Z.Y. Sapiro 
|16j for a detailed presentation of these generalized spherical harmonics. 
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There exists thus a family of function with the indexes (I, n) running in the set I = {(£, n), Z — |g| G 
N, I - \n\ G N} which is a Hilbert basis of L 2 (S 2 ; C 2 ) and such that V(Z, n) G X, 

The Hilbert space H. can then be decomposed into the infinite direct sum 

H= [L 2 (R X ;C 2 )®^] := H ln , 

(Z,n)eZ (J,n)£I 

where H/„ = L 2 (E K ;C 2 ) ® is identified with L 2 (R;C 2 ). We obtain the orthogonal decomposition for 
the Hamiltonian H 

#=0 H ln , 

(Ln)el 

with 

H ln :=H lHln = T 1 D x +a l (x)T 2 , (2.11) 

and ai(x) = —a(x)(l + |). The operator iJ'" is a selfadjoint operator on Tti n with domain D{H ln ) = 
i? 1 (IR;C 2 ). Note also that we use the following representation for the Dirac matrices T 1 and T 2 . 

In this paper it will be enough to restrict our study to a fixed harmonics. To simplify notation we 
thus simply write H, H and a{x) instead of Hi n , H ln and a/(x) respectively. 

We summarize now the direct scattering results obtained in [231 HH E|- We define H = r 1 !)^ as the 
free Hamiltonian. It is clearly a selfadjoint operator on TC with domain D(H n ) = i? 1 (K;C 2 ). The main 
results concerning the spectral properties of the pair (H, Ho) are the following. The Hamiltonian H has 
no eigenvalue and no singular continuous spectrum, i.e. 

a pp (H)=<D, a sc (H)=®, (2.13) 

and the standard wave operators 

W ± = s- lim e itH e- itHo , (2.14) 

t — '±OC 

exist and are complete on TL. Let us make a few comments on these results. It follows directly from 
(|2.13|) and (|2. 14[1 that H has purely absolutely continuous spectrum and thus all states scatter in the 
asymptotic regions of spacetime where they obey a simpler equation given by the comparison dynamics 
e ltH ° . Notice that this comparison dynamics is particularly simple for massless Dirac fields. Due to the 
diagonal form of the matrix r , it becomes simply a system of transport equations with velocity +1 for 
the second component of the spinor and velocity —1 for the first. At last, observe that the potential a 
has the following asymptotics as x — > ±oo: 

3k >0, |a(ar)| = 0(e KX ), x -» -oo, 

\a(x)\ = 0(|x| _1 ), x +oo. 

Eventhough a is apparently a long-range potential having Coulomb decay at +00, we don't need to modify 
the comparison dynamics in the definition of the wave operators. We refer to [7], Section 7 for a proof of 
this point. 
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The scattering operator S is defined by 

S=(W + )*W~. (2.15) 

We can simplify the operator S by considering incoming and outgoing initial data separateljU. In order 
to define these subspaces of initial data precisely, let us introduce the asymptotic velocity operators 

P± = s-C oa - lim e ttH o^e- itH °, 

and 

P ± = s-C oa - lim e ltH -e' ltH . 

t^±OC t 

We refer to [5] for a detailed presentation of these operators and their usefulness in scattering theory. By 
a direct calculation, we get easily P^ = T 1 while it was shown in [7] that the operators P ± exist, have 
also spectra equal to {—1, +1} and satisfy the intertwining relations 

P ± W ± = W ± T 1 . (2.16) 

Let us denote P ln = l R -(r 1 ) and P ou tlK+(r 1 ) the projections onto the negative and positive spectral 
subspaces of the "free" asymptotic velocity T 1 . Then we define Hf n — Pi a H and H® ut = P ou _tH and refer 
to these spaces as incoming and outgoing initial data for the free dynamicfjf). In a similar way, we define 
H in = 1 K -(P ± )H and H^ ut — 1 K +(P ± )H the spaces of incoming and outgoing initial data for the full 
dynamics when t — > ±oo. In our particular case it turns out that we have equalities between some of 
these spaces. Precisely we have 

Lemma 2.1 With the same notations as above, we have 

7~tout = W ou( =: Ti-out- 

Proof: We only prove the case "out" since the proof for "in" is similar. We shall use the following 
characterization of the asymptotic velocity operators obtained in [7] 

P ± = s-C oc ~ lim e itH T l e- itH . (2.17) 

t— >-±oo 

We compute l R _(p-)l R+ (P+). By (|2~T7) we get 



1 R - (P-)1 R+ (P+) = 8 - A lim e~ ltH 1 R _ (r^e 2 "^ 1 R+ (r 1 



^e~ UH . 



Since T 1 !! = HT 1 with H = T 1 ^ - a(x)T 2 , we have 



1 K -(P-)1 K+ (P+) = a- lim e-^e^lR-CT^lB+^Je- 

= o. 



2 The terms incoming and outgoing data have to be understood here in reference to the black hole, i.e. as the data whose 
associated solutions propagate, when t increases, in direction of the black hole (incoming data) or escape from it (outgoing 
data) 

3 Observe that Hf n (resp. H° at ) are simply the spinors with nonvanishing first component (resp. second component) 
only 



(i 



Hence 1 R +(P + ) = (l K - (P ) + l E +(P ))l R +(P + )l R +(P )1 B + (P + ) from which we deduce immediately 
that Ti.^ nt C Tto Ut - Analogously we can prove H~ ut C TQ ui and thus the equality TQ ut — TC~ ut holds. 

<> 

Thanks to the intertwining relations (|2.16p and Lemma \2. II we see that the wave operators W ± are 
partial isometries from H? n ^ out into 7ii n / ou t- We define 

W t t = s- lim e ttH e- ttHo P ouU 

£— >±oo 

W± = s - lim e ltH e- ltH " P in , 

t — >±oo 

the corresponding outgoing and incoming wave operators. Remark here that by the relations £/o1r± (T 1 ) = 
±D a; l R ±(r 1 ) the Hamiltonian Hq becomes "scalar" when projected on W? n / out . Hence we can simplify 
the expression of the outgoing and incoming wave operators by 

W± t = s- lim e itH e- itD *P out , 
W± =a- lim e ltH e ltD * P in . 

t — >±00 

This minor remark will turn out to be important in the next section. Now from the definition (|2. 15|) of 
the scattering operator, we deduce that S leaves invariant the incoming and outg oing subspaces 7~^i n / ou t 
and thus can be decomposed into the direct sum 

S = S in ®S out , (2.18) 

where 

S in = (W+)*WT, S out = (W + ut yW ~ t . (2.19) 
Clearly S- ln and S Q ut ar e isometries on 7i? n and H® ut . 



3 The inverse problem at high energy 
3.1 The main result and the strategy of the proof 

We are now ready to study the inverse problem for massless Dirac fields satisfying (|2.8|) . The question 
we adress is: can we determine the mass M and the charge Q of the black hole from the knowledge of 
the scattering operator S? The main result of this paper is summarized in the following Theorem 

Theorem 3.1 Assume that the scattering matrix S ou t (or Si n ) is known. Then the corresponding poten- 
tial a(x) can be entirely determined. Moreover the mass M and the square of the charge Q 2 of the black 
hole can be recovered by the formulas 

x — x 3 a 2 (x) .„ „ , 

M= lim —U., (3.1) 

x — *+oo z 

Q 2 = lim (x 4 a 2 (x) - x 2 + 2Mx). (3.2) 

x — >+oo 
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The strategy we adopt to prove this Theorem is based on a high-energy asymptotic expansion of S ou t , 
a well-known technique initially developed in the case of multidimensional Schrodinger operators by Enss, 
Weder [IT]. This method can be used to study Hamiltonians with electric and magnetic potentials [JJ, 
the Dirac equation [21] and Stark Hamiltonians [3TJ [27] . 

Precisely, we consider ipi,ijj2 £ ?^out such, that ipi,ip2 & Co°(lR; C 4 ) and we define the map 

Font (A) =< Soute^Vi, e lA ^ 2 >, (3.3) 

for A G K. Following the ideas of 25, 26 , we shall obtain an asymptotic expansion of S(X) when A — > oo 
with the exact expression of the first three leading terms. Exactly we get the following reconstruction 
formula: 

Theorem 3.2 (Reconstruction formula) Let tpi, ip 2 € TiPout su °h that ipi,ip2 & Cq°(R; C 2 ). Then for 
A large, we obtain 

F out (X) =< >+j< il>\,L x {x,D x )il)2 > < ifJi,L 2 (x,D x )ip 2 > + 0(\- 3 ), (3.4) 

where Lj(x,D x ) are differential operators given by 

L 1 {x,D x )=L 1 = 2 > , 
2r + 



8r± 2r 4 



and r+ = M+ M 2 - Q 2 



Let us admit temporarily this result and use it to prove the main theorem. If we assume that the 
scattering matrix 5 ou t is known, then we can use inductively the high-energy expansion (|3.4j) to recover 
first, the radius r+ of the exterior event horizon (term of order A -1 ) and second, the potential a(x). 
Indeed, the term of order A -2 (since r+ is already determined) gives the quantity 

<^Vi,V>2>. (3.5) 

Since Q3.5P can be determined for any if>i,ip2 in a dense subset in Tt, we recover completely the potential 
ai(x) and thus a(x). Now the formulae for the mass p.l|) and for the charge (|3.2p follow directly from 
the definition of a(x) 

a 2 {x) = ^±21, 

and the fact that 

r(x) ~ x , x — > +oo 



an immediate consequence of (|2.5[) . This ends the proof of our main Theorem and proves that observers 
static at infinity who can measure the scattering matrix are able to recover the mass of the black hole 
and its charge, up to a sign. Note that the first term of the expansion (|3.4j) permits us to determine the 
radius r + of the exterior event horizon, an information of physical relevance. 
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The main step in the proof of Theorem 13. II is clearly the reconstruction formula [33] stated above. In 
order to prove it, we first rewrite F out (A) as follows 

F out (X) = <S oxlt e iXx ^ u e iXx ^ 2 >, 

= <W- t e lX ^ 1: W+ t e lXx ^>, 

= <W- ut (X)^uW+ t (\)^ 2 >, (3.6) 



where 
and 



W±(X) = s- lim e *ffW e -«(o.+A)p a 

t— »±oo 



H(X) =T\D x + X) + a{x)T 2 . 
In order to obtain an asymptotic expansion of the operators W out (A), we follow the procedure exposed 



in [25] [26] , procedure inspired by the well-known Isozaki-Kitada method 20J developed in the setting of 
long-range stationary scattering theory. It consists simply in replacing the wave operators W oui (X) by 
"well-chosen" energy modifiers J ( ^ lt (A), defined as Fourier Integral Operators (FIO) with explicit phase 
and amplitude. Well-chosen here means practically that we look for J^ nt (X) satisfying for |A| large 

W£ t (AW= lim e « ff W4(A)e-«( D ' + ^, (3.7) 

t— >±oo 

and 

||(W* t (A) - Jt^mW = 0(X- 3 ), (3.8) 

for a fixed ?/> S T^out such that $ e Co°(IR; C 2 ). In particular if we manage to construct J^t(^) satisfying 
(pT8]) then we obtain by f[3T6|) 

F out (A) = < J- ut (A)^, J+ t (A)V>2 > + 0(A- 3 ), 

from which we can calculate the first terms of the asymptotics easily. 

The rest of this paper is devoted to the construction of the modifiers J^ t (A) and to the proof of 
Theorem 13.21 For simplicity of notations we omit in the next subsections the underscript "out" when 
naming the wave operators W^ t and modifiers J out . 



3.2 Construction of the modifier 

Let us give a hint on how to construct the modifiers J (A) a priori defined as an FIO with "scalar" phase 
l ^ ± {x 1 £, A) and "matrix- valued" amplitude p^(x, £, A). If we assume that (|3.7[) is true then we can write 

/■±oo 

(W ± (\)-J ± (\))Tp = i / e ttH ^C ± (X)e- lt(D - +X) ^dt, 
Jo 

where 

C±(A) := H(X)J ± (X) - J^AXAr + A), (3.9) 
is an FIO with phase (p^(x, £, A) and amplitude c ± (a;, £, A). Hence we get the simple estimate 

/>±oo 

\\(W ± (X)-J ± (Xm< IIC+CA) e- itD ^\\dt. (3.10) 
Jo 

In order that (|3.8|) be true it is then clear that the FIO C(A) has to be "small" in some sense. Precisely 
we shall need that the amplitude cF^x, £, A) be short-range in the variable x and of order A -3 in the 
variable A. 
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3.2.1 Construction of the modifiers at fixed energy 

We first look at the problem at fixed energy (i.e. we take A = in the previous formulae). Hence we aim 
to construct modifiers with scalar phase <p (x, £) and matrix-valued amplitude p^ix, £) such that the 
amplitude c ± (a;, £) of the operator C ± = HJ^ — J^D X be short-range in x. We follow here the treatment 
given by Gatel and Yafaev in [15] where a similar problem for massive Dirac Hamiltonians in Minkowski 
spacetime was considered. 

The operator is clearly an FIO with phase ip ± (x,£ l ) and amplitude 

c^z.o = b ± (x,Op ± (x,o - ir^pifoo, (3.ii) 

where 

B ± (x,0 = T 1 d x i P ± (x,0 + a(x)T 2 - ^. (3.12) 

As usual we look for a phase ^> ± close to x£ and an amplitude p ± close to 1. So the term d x p ± in (13 . 1 1 [> 
should be short-range et can be neglected. Moreover if p^ was exactly 1 we would have to determine (p^ 
so that 

B ± =0. (3.13) 

A direct calculation of S ± = leads to a matrix-valued phase whereas we look for a scalar one. But if 
(p^ 1 satisfies (|3.13p then it should also satisfy the equation B 2 = 0. Using the anticommutation properties 
of the Dirac matrices we have 

{B±f{x, = {d x ^f + a 2 (x) -e- 2^B ± (x, £) = 0. (3.14) 

Using (|3.13p again, we eventually get the scalar equation 

f^foO :={d 3lV > ± ) 2 +a\x)-^ = 0, (3.15) 

in fact an eikonal equation of Schrodinger type. Let us set ip — x^ + (fr^ix,^) where ^(x,^) 

should be a priori small in the variable x. Then we must solve 2^d x <j) ± + (9 2; ± ) 2 + a 2 (x) = 0. Neglecting 
(d x 4r^) 2 in this last equation, we obtain 

2£0 !B ± +a 2 (x) = 0. (3.16) 
For £ ^ 0, we get the following two solutions of (|3.16[) by the formulae 

1 f ±oa 

^( x ^) = -l a 2 (x + s)ds, (3.17) 



o 



and with this choice, we obtain for £ =/= 0, 

r ± (x,0 = ^a*(x). (3.18) 

In our derivation of the phase, it is important to keep in mind that we didn't find an approximate 
solution of (|3.13p but instead of (|3.14p . Therefore we cannot expect to take p ± = 1 as a first approximation 
and we have to work a bit more. In order to construct the amplitude p ± we follow here particularly closely 
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So we look for p^ such that B ± p ± is as small as possible. What we know is that (B^) 2 is already 
small thanks to our choice of tp^. Hence let's try to find a relation between and (B^) 2 . Observe first 
that B^ have the following expression on 

B ± {x, = B(x, = B (x, + 2£,K{x 7 0, (3.19) 

where 

B (x,0 = -2£P ln , (3.20) 

K(x,0 = ±(-^Ti+a(x)T 2 ). (3.21) 

In particular, B^ don't depend on ± in T ± . Taking the square of (|3 . 19[) we get 

B 2 = B 2 + 2£B K + 2£KB (3.22) 

But from (|3~20|) we see that Bq = -2£B . Whence (|3~22|) becomes 

B 2 = -2£B (1- K)+2£KB. (3.23) 

If we add 2£B to both sides of the equality in (|3.23|) then we obtain by (|3.14|) . 

r[x, := B 2 + 2£B = -2£B (l - K) + 2£(1 + K)B, (3.24) 

or equivalently, if we isolate B on the left hand side 

2£(1 + K)B = r(x,0 - 4£ 2 P in (l - K). (3.25) 

It follows immediately from the definition oiK that (l±K) are inversible for I := { ( e I : | £ 1^ }> 
R 3> 1. In consequence we can write for 

B{1 - K)- 1 = -Ul + K)- X r{x, 0(1 - K)- 1 - 2£P in . (3.26) 

The first term in the right hand side of (|3.26[) is small thanks to our choice of phase but the second one 
is not. We choose p^ in such a way that it cancels this term. According to (|3.26p . a natural choice for 
p ± is thus 

p ± (x,C) =p(x,0 = (1 - *0 _1 ^out, (3.27) 

for which we have 

q(x,0 :=B(x,t)p(x,0=^(l + K)- 1 r(x,Z)(l-K)- 1 . (3.28) 
Note that p defined by l|3.27p is independent of ±. 

Let us summarize the situation at this stage. For£ ^ 0, we have defined the phase ^(x, £) =x^+0 ± (x,O 
by (|3.17|) and for £ e X, the amplitude p is given by (|3.27|) . Directly from the definitions, the following 
estimates hold. 

Lemma 3.1 (Estimates on the phase and the amplitude) 

V(x,OeRxX,Va,peN, |c£ctV(a;,£)| < C afi (1+ | x \ + \£ |)( 2 "«-« + . (3.29) 
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V(x,£)eKxX, \dl^(x,0-xO\<-^ (3.30) 



V(a;,OeR ± xX,Va,/3eN,, ^^(x, 01 < C a p (x)- 1 -"®- 1 -". (3.31) 

V(i,0 elx X,Va,/3 G N, |c£cf K(x,OI < C Q/3 (x)- 1 -^ ^ (3.32) 

V(i,0 elx X,Va, (3 G N, |d£9f (p(x,0 - P out )| < C Q/3 (x)- 1 - Q (e)- 1 ^' 3 . (3.33) 

V(i,OeixI>,/3€N, |0°afr(a;,OI < C Q/3 (aT 4 ~ a (O^- (3.34) 

V(x,£) GExX,Va,/3GN, |5^g(x,£)| < ^ (x)- 4 ""^)- 3 ^. (3.35) 

Thanks to (|3.29| . (|3.30p and (|3.33p . for R large enough, we can define precisely our modifier J ± as a 
bounded operator on H ([29], Corollary IV. 22): let 9 G C*°°(R) be a cutoff function such that 0(£) = 
if |£| < | and 0(£) = 1 if |£| > 1. For R large enough, is the Fourier Integral Operator with phase 
9? ± (x,^) and amplitude 

P(x,0=p(x,0 9(j i ) (3.36) 

Moreover, we have 

Proposition 3.1 For any -0 G T~P out such that -0 € C^°(X;C 2 ), we ftawe 

W ± i)= lim e <tff J ± e -<t2J -^. (3.37) 

t— ►ioo 

Proof: We only consider the case (+). It is clear that it suffices to show 

f^ + oo V R J 

Let ?7 G C°°(R) such that rj(x) = 1 if x > 1 and = if x < 0. We write e~ ltD *^ as : 
Clearly, (1 - <q{x)) e~ ltD *V = e~ ltD * (l-<q{x + t)) tp -> on H. 

Now, using (|3.29[) . (|3.3ip , (|3.33|) and the standard pseudodiffcrential calculus, we see that 
T := ( J + — 9(—^-)P out ) rj{x) is a ^Ao with symbol fc + (x,£) which satisfies 



(3.38) 



R 

V(x,£) G R x R,Va,/3 G N, |d«^fc + (x,£)| < C Q/3 (x)- 1 - ^) -1- ^. (3-40) 



Then T is a compact operator on TL and since e UD ^ip — > weakly, Proposition ^. H is proved. 
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3.2.2 Construction of the modifiers at high energy 

Let's come back to the construction of ^ ± (A) the modifiers at high energy. Recall that we look for 
modifiers that satisfy (13. 7\i and (|3.8p . Proposition 13.11 suggests to construct J ± (A) close to e~ l X J e . 

Clearly, e~ lXx J ± e lXx are the FIOs with phases ^{x, £, A) = x£+ ^ 2 (g+A) Jo ±0 ° a2 ( x + s )ds\ and amplitude 

P(a:,£ + A). Since the modifiers J ± (A) will be applied on functions tp having compact support in the 
Fourier variable, the cut-off in the £ variable in the definition of P(x,£ + A) disappears for A sufficiently 
large enough. 

Thus, if we take exactly J ± (A) = e~ lXx J^e %Xx as the modifiers, according to p. lip , for £ in a compact 
set and A> 1, the amplitude of the operator C ± (A) becomes 

<*{x, f , A) = c(x, £, A) = P(x, £ + A)p(x, £ + A) — ir^pfo £ + A). 

We want the amplitude c(x, £, A) to be of order 0((a;)~ 2 A~ 3 ). By (|3.35p the first term P(x, £+A)p(x, £+A) 
has order 0((a;)~ 4 A~ 3 ). By (|3.33|) the second term d x p(x,£ + A) has order 0((a;)~ 2 A _1 ) which is not a 
sufficient decay in A for our purpose. 

So, we have to work a bit more. We look for the modifiers J ± (A), with A large enough, as FIOs with 
phases y ± (x, £, A) and with a new amplitude P(x, £, A) under the form 

P(x,S,X) = (p(x,S + *) + ^(Ptaki(x,S)+p(x,Z + \)h(x,zj) +^P m fc 2 (i,0) 9(0, 

where g € Cq°(M) with g = 1 on Supp ip, the functions k\,k2,h (that can be matrix-valued) will be 
short-range of order 0((x)~ 2 ). Observe that there are two different types of correcting terms in this new 
amplitude: on one hand the terms which contain p(x, £ + A) ~ P ou t (when A — > oo) are correcting terms 
that "live" in H.® ut ; on the other hand the terms which contain P; n are clearly correcting terms that "live" 
in . Now, noting that 

B(x, £ + A) = -2(£ + A)P in + M(x, £ + A), M(x, £ + A) - - ^M- T 1 + a(a:)r 2 , 

the amplitude c(x, £, A) can be written for £ 6 Supp ip as 

c(x,£,A) = P(x,£ + X)P(x,£ + A) - %T 1 d x P{x,C + A), 
= B(x,Z + X)p{x,Z + X)-ir 1 d x p(x,£ + X) 

+ ^ ( - 2(£ + A)P in fc!(x, + M (as, £ + A)P in /ci(x, £) - il^ftAfcifo £)) (3.41) 
+ ^ (p(x, £ + A)p(x, £ + A)/i (x, - iT 1 (d x p(x, £ + A)/i + p(x, £ + X)d x h {x, £))) 
+^3 ( - 2(£ + A)P in fc 2 (x, + M(x, £ + X)P in k 2 (x, £) - zr 1 P in a,fc 2 (x, o) ■ 
Directly from the definitions we get the following asymptotics when A — > oo 

d xP (x,z + x) = +^r 2 p out -^Mr 2 p out + o((x)- 2 A- 3 ), 

p(a:,£ + A) = Pout+O^^A- 1 ), 
P(x,£ + A);p(x,£ + A) = 0((a;)- 4 A- 3 ). 
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This leads to 

c(x, £, A) = -i^T'T 2 P oat + i^W T i r2pout _ 2 _ 2| p;nfci 

+ ^r 2 P in fci - i^^Pi^fei - i^PoMi - ^P in k 2 (3.42) 
+0((x)- 2 A- 3 ). 

But the anticommutation properties of the Dirac matrices entail the relations r 1 r 2 P out = — p n r 2 and 
r 2 P in = P out r 2 . Hence the amplitude c(,£,A) takes the form 

1 D (; a '( X )-n2 



c(x,t,\) = -P in (i^r 2 -2fc!) (3.43) 



r 2 - 2^fci + id x k x - 2fc 2 ) + ^P out (a(x)T 2 ki - id x h 



A 2 V " 2 * " s '^ ' ' A 2 ' 

+0((a;)- 2 A- 3 ). 

Using (I3.43| we can cancel the terms of order less than 0(A~ 3 ) inductively. We choose 

k 1 (x,Z)=i^I a , (3.44) 
which is obviously short-range of order 0{(x)~ 2 ) and cancel the first term in (|3.43p . It remains then 

c(a5,C,A) = TjPn^-« — 2 — r ~* — 4 — r 4 — ~ 2fc2 J ' 3 ' 45 ) 



1 n /.a'(a;)a(x) ,„ , \ ._o,_°, 
^ -Pout ^ 4 ~ *$r*i j + 0«a) 2 A 3 ), 



for £ S Supp if). We choose 



and 



h(x,0 = ^P-, (3-46) 

k 2 (x,o = -i^-tr 2 -—^r 2 , (3.47) 

The three correcting terms k\,li, ki are short-range of order 0((x)~ 2 ) for £ in a compact set and cancel 
the two first terms in (|3.45|) . This leads to 

V(x,£)eR 2 ,Va,/?€N, \d^c(x,Z,\)\ < C a0 (x)~ 2 - a A~ 3 . (3.48) 

Eventually the new amplitude P(x,^, A) of J (A) is defined as 

P(x, f, A) = (p(z, € + A) + 1 (i^P in r 2 + £ + A)^ 

+ ^(-^P^-^p ln r 2 )) 

or equivalently 

f, A) = (p(x, £ + A) + 1 (^r 2 F out + £ + A)^M 
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1 I »)t, 2d 



A , , • , r^P out - -^P out j j (3.49) 
Mimicking the proof of Proposition 3.1, we have 

Lemma 3.2 For any ip € H^ut such that € C^°(R) and for X large, we have 

W ± (X)ip = lim e Jtff ( A )j ± (A)e-' it ^ +A V. (3.50) 

i — »±oo 

Now, we can state the main result of this section: 
Lemma 3.3 When A tends to infinity, the following estimate holds: 

\\(W ± (X)-J ± (\m=0(X- s ). 
Proof: Let us consider the case (+). Thanks to (|3.10p . we have to estimate ||D + (i, X)ip\\ where D + (t, A) 



is the FIO with phase 



1 



(p + (x,£,X,t) = x£+ 2 (^ + - x) J a 2 (s) ds 



+ OC 



(which is uniformly bounded with respect to t) and with amplitude c + (x + t, £, A). 

Let x € C§°(R) be a cut-off function defined by x(x) = 1 if | x |< §, x(%) = if | x |> 1 and let C = l~X- 
We have 

2t 2r 
A)V|| < ||x( T ) J D + (t, AMI + ||C( T )^ + (^ A)^|| := (1) + (2). (3.51) 

First, we estimate the contribution (1). On Supp x(^r)i I ^ + £ |> |. Using (|3.48[) and the continuity of 
FIOs, we have 

2r 

||X( T )I> + (*,A)|| = 0(<i>- 2 A" 3 ). (3.52) 

Now, let us estimate the contribution (2). It is clear that 

2r Rr 2r 8r 

(2) < || C (_)D+(t,A)C( T )VII + \\C( T )D + (t,X)x( T M := (a) + (6). (3.53) 

On Supp C(^), | x |> j s . Since ^ S 5(R), HCCfOV'll = 0(< < >~ W ) for all JV > 0. Thus, using the 
continuity of FIOs again, we obtain (a) = 0(< t >~ N A -3 ). 

Now, we evaluate the contribution (6) by using a standard non-stationary phase argument. We have 



C&D+(t,X)x(^Mx) = (2tt)- 2 C(^) / e^ + ^-y^c(x + t,t\)x( 8 -^)ip(y)d{;dy, (3.54) 

t I t ,/ B 2 / 

where 



$ + (t,s,frg,A) = (a:-g)£+ ' / a 2 (s) ds. (3.55) 



1 

In order to investigate possible critical point, we calculate 



dt* + (t,x,y,Z,\)=x-y- I a 2 {s)ds. (3.56) 

AS + A J Jx+t 

On Supp C(x)' I x I— I ana - 011 ^"PP x{.~t)-> I 2/ I— f • Moreover the integral which appears on the (RHS) 
of (|3.56p is uniformly bounded with respect to X,x and t. Thus, | d^ + (t,x,y,£ y ,X) |> C(l + t+ \ x |). 
We conclude by a standard argument of non-stationary phase that for all N > 0, (b) = O (< t >~ N A~ 3 ). 
Integrating over (0, +oo), we obtain Lemma [3~3l 
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3.3 High energy asymptotics of the scattering operator 

In this section we use the previous construction of the modifiers (A) to prove the reconstruction formula 
(|3.4[) . Recall that we want to find the asymptotics as A — ► oo of the function (see (|3.6[) ) 



F{\) =< W-(A)Vi, W + {\)^ 2 >, 

where have compact support. Using Lemma l3.3i we get 

F(\) =< J-(A)V>i, J + (A)^ 2 > + 0(\- 3 ). 

Let us first give the asymptotic expansion of J ± (A) at high energy. By the constructions above, we have 
for A large 

J±(A)^ = ^| R ^ ± (^)( p(x ,e + A) + ^(z^r 2 +P (x,e + A)^M) (3.57) 

Thus J ± (A) can be seen as pseudodifferential operators whose symbols j -± are given by 

3 (x,£,X) = e 2<«+A)Jo \p[x,£ + A) + —\i — - — T +p(x,^ + A) — - — ) 

1 / £a'(x)„<> a"(x)„o\\ 
+^(-^ r "4 ir ))- 

Hence, using Taylor expansion of e at t = and using the expansion 

+ A) = (i + ^r 2 - ^r 2 )p out + o(\-% 

we get 

(l + ^ r3 _^ r2 + ! ^, r2 + ^ +0(A _ 3)) , 

where we denoted -/^(x) = Jq ±0 ° a 2 (^ + s)ds. If we put together the terms with same order, we obtain 

3 ( a> e,A) = i + (_ + _r j + + 

We denote i? the differential operators 

i / il ± D x (Z*) 2 a 2 (a;)\ (iaix)!^ a(x)D x a'(x)\ 9 ,„ „ . 

fl = ( — - s + -F) + (4 4^ + ^) r < (3 - 58) 

as a shorthand. Then the high energy asymptotics of J (A) are 

j± m = $ + 1 + ^r 2 )v. + + 0(A- 3 ). (3.59) 
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We now can prove Theorem 13. 21 We have 
F(X) = < J"(A)^i, J+(A)^ 2 > +0(\- 3 ) 

4( < > + < *.« + * > ) + <(£ + 1? r2 K(£ + tM*}' 

M JV 
= < ^>i,Va > +— + 

Let's compute the first order term M. Since r 2, !/>i G 7Y" n , the following terms vanish 

< rVi, Va >=< V>i, r 2 V>2 >= 0. (3.60) 

In consequence 

M = < i/> u Ul+ - r)i/>a > = % - I a 2 (s)ds < ^i,V2 > ■ 
But a direct calculation using (|2.2| and (|2.4|) shows that / a 2 (s)ds = — . Whence 

M= \ 2> <Vi,^2>- 3.61 
2r+ 

We compute now the second order term N. Notice first that by (|3.58p and (|3.60p . we get 

<R^ 2 >=({-— Sr"^)^ 2 /' 

and 

/fil~ a(x) o\ , /*^ + a ( a; )^2\ , \ , l T -r+, , a2 i x ) i 

{{— + JV*, (— + -^r 2 ]^) =< -/ /+V-2 > + < Vi, -j^V* > • 

Hence we obtain 

at = - < ^, - j-)a^2 > +| < Vi, (i + - n 2 ^ > + <ih, ^V> 2 > • 
f ( l + 1 ) 2 

Using / a 2 (s)ds = — again, 

Jr r + 

N = - l l±j)l < ^,D X ^ 2 > + < Vi, 4^2 > + ^#^ < Mi > ■ (3-62) 
2r+ 2 8r+ 

The two leading terms (I3.61|) and (|3.62|) give exactly the reconstruction formula (|3.4p . Hence our main 
result is proved. 
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